Noise Effects on the Complex Patterns of Abnormal Heartbeats 
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Patients at high risk for sudden death often exhibit complex heart rhythms in which abnormal 
heartbeats are interspersed with normal heartbeats. We analyze such a complex rhythm in a single 
patient over a 12-hour period and show that the rhythm can be described by a theoretical model 
consisting of two interacting oscillators with stochastic elements. By varying the magnitude of the 
noise, we show that for an intermediate level of noise, the model gives best agreement with key 
statistical features of the dynamics. 

PACS numbers: 5.40. Ca, 5.45. Tp, 87.19.Hh, 89.75.Kd 



Individuals with frequent abnormal heartbeats (Fig. [l]) 
may be at high risk for sudden cardiac death JO . Such ab- 
normal heart rhythms often have a random appearance. 
Attempts have been made to analyze these rhythms by 
inspecting short data strips @|| and matching them beat 
by beat to various models |^|]]. Other approaches 
characterize statistical properties in longer records of up 
to several hours. However, the mechanisms underlying 
these abnormal rhythms and their changes over long pe- 
riods of time remain elusive. Here we show that a the- 
oretical model consisting of two coupled oscillators J^-fJ] 
describes the observed patterns of abnormal heartbeats 
in one clinical case provided we introduce noise to the 
periods of the oscillators. This approach may generalize 
to the analysis of the underlying mechanism of a large 
number of records with complex patterns of abnormal 
heartbeats. 

We consider a continuous 12-hour segment of the 
ambulatory electrocardiographic record of an individ- 
ual with heart failure and frequent abnormal heartbeats 
(Fig. [l]). The normal sinus heartbeats, S-beats, arise from 
activity in the sinus node, the normal pacemaker of the 
heart. The time intervals between the normal beats ap- 
pear to be periodic in this short tracing, but they do 
fluctuate during the 12-hour period. The abnormal ven- 
tricular beats, V-beats, arise in the lower chambers of 
the heart, the ventricles. Although the timing between 
the V-beats appears to be irregular, the histogram of 
the interventricular time intervals, the time intervals be- 
tween consecutive V-beats, consists of equidistant peaks 
implying that the interventricular time intervals are mul- 
tiples of a fixed number (Fig. [I]). This result is consistent 
with the possibility that the heart rhythm arises from a 
competition between two oscillators: the normal sinus os- 
cillator, and an abnormal ventricular oscillator with pe- 
riods T$ and Ty, respectively. Such rhythms are called 
parasystolic. 

After each S-beat or V-beat, there is a time period, 



called the refractory time 9, during which all other heart- 
beats are blocked. Thus, all V-beats must occur at least 
a time 9 after an S-beat. Because in our case 9/T$ > 1/2, 
such a V-beat will block the following S-beat |l(]] . 
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FIG. 1. Upper panel: Electrocardiogram over 6.6 s of a pa- 
tient with heart failure. Sinus (S) and ventricular (V) beats 
differ in morphology. S-beats are separated by an interval 
Ts- The time interval between the S-beats with an inter- 
vening V-beat is approximately 2 Ts implying that an S-beat 
has been blocked (downward white arrow). The interventric- 
ular time interval Ty, determined from the histogram in the 
lower panel, is also indicated and is consistent with a blocked 
V-beat (upward white arrow). Lower panel: The histogram 
of the time intervals between consecutive abnormal V-beats 
for the 12-hour record from which the segment in the upper 
panel was extracted (« 80 000 beats with « 5000 V-beats). 
The time interval between the peaks is 1.67 s (« 3Ts). The 
inset is an enlargement of the peak at 5 s. The sampling rate 
of the underlying electrocardiographic time series is 128 Hz. 

To describe the parasystole mechanism we introduce the 
phase <f>i of the ith ventricular beat in the sinus cycle, i.e., 
the time interval between the zth ventricular beat and the 
previous sinus beat normalized by Ts [3™]- Successive 
values of fa are determined by iterating the difference 
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equation: 
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where a ventricular beat is expressed if 4>i+i > @/Ts- 
For a fixed irrational ratio Ty /Tg (incommensurate pe- 
riods), the phase of the V-oscillator with respect to the 
S-oscillator will, for sufficiently long times, be equally dis- 
tributed in the interval [0,1]. Interestingly, the sequence 
of the number of intervening S-beats between consecutive 
V-beats, called the NIB sequence, consists of only three 
different values (NIB-triplet). The values in these NIB- 
triplets change with Ty /Ts and 9/Ts |0. Examples are 
shown in Fig. |^a. 

The uniform distribution of the phases of the V-beats 
allows us to derive the fraction of V-beats. Call ny, the 
number of V-beats observed on the electrocardiogram in 
a time interval T, and N the total number of S-beats and 
V-beats during the same time interval. The total number 
of observed and blocked V-beats in T is T /Ty and the 
total number of observed beats is N — T/Tg. The frac- 
tion of observed V-beats, i.e. those which occur outside 
of the refractory period is (Ts — 6) /Ts ■ It follows that 
ny = (T s - 0)/T s x T/T v , and consequently 



ny/N = (T s -9)/T v . 



(2) 



Perturbation of the timing of the V-beats by Gaussian 
noise does not affect this distribution, because the uni- 
form distribution of the phases of the V-beats remains 
unchanged. 

To illustrate the properties of the unperturbed model, 
we iterate Eq. (|l|) using typical parameter values for the 
patient under consideration. Histograms showing the 
NIB values for three sets of parameters are shown in Fig. 
|a. 

We now compare the results of Eq. ([!]) with the clin- 
ical data. Since the value of Ts fluctuates during the 
record, we combine the data from different times of 
the day during which Ts falls in a fixed 10 ms range. 
Fig. ||b shows the histograms of NIB values for three 
different values of T s . For T 5 =0.61 s (Ty /T s < 3) 
and Tg=0.51 s (Ty /Ts > 3), we find the same NIB- 
triplets as in the deterministic model (Fig. ||a). How- 
ever, additional peaks in the data (Fig. ||b) contradict the 
"rule" of only three NIB-values in the purely determin- 
istic model. Furthermore, for T s = 0.55 s (T v /T s « 3), 
fewer beats with NIB=2 are present in the data, and new 
NIB-values appear largely corresponding to the sequence 
5, 8, 11 . . . = 'in — 1, where n > 2 is an integer p2[ . 

There are also significant discrepancies between the 
predicted fraction of V-beats from Eq. (g) and the ob- 
served values. Fig. [|a shows that, for Ty /Ts ~ 3, the 
model of Eq. ([l]) predicts many more V-beats than are 
found. 
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FIG. 2. Generation of NIB-triplets for d/T s = 0.5 and 
two different values of Tv/Ts ■ The boxes are the S-beats 
including their refractory times. An empty box is a blocked 
beat, and a dashed box an expressed beat. The expressed 
V-beats are symbolized by black bars, the blocked ones by 
empty bars. The S-beat after each expressed abnormal beat 
is blocked since 9/Ts > 0.5. (a) Deterministic model with- 
out coupling. In the upper panel, the ratio of the peri- 
ods Tv/Ts ~ 3. The phase of the V-beat with respect to 
the S-beats is therefore changing slowly, and the NIB is 2 
for several periods until the V-beat falls into the refractory 
time of the S-beats and is blocked for several periods. For 
Tv/Ts « 3.3 the NIB-triplet 2,6,9 is generated, (b) Model 
with coupling and noise. Coupling shifts the V-beats into the 
refractory time for Tv/Ts ~ 3. For Tv/Ts ~ 3.3 the coupling 
has minimal effects. The addition of noise (arrows) changes 
the NIB-sequences by moving some V-beats into or out of 
the refractory time. For Tv/Ts ~ 3, this combines two con- 
secutive appearances of NIB=2 into NIB=5 (8,11,. . .). For 
Tv/Ts — 3.3, combinations of NIB=2 with one of the other 
two values, 6 or 9, generates NIB=9 or 12. 

We propose that the pattern of V-beats in the data 
may be understood by assuming that the timings of the 
V-beats are not strictly periodic as in Eq. ([j]), but that 
they are described by a stochastic difference equation |ll| 

, T v . f(4>i,Ts,T v ) 

+ - + 



T s 



modi, (3) 



where 77 is a Gaussian random variable distributed 
around 0, and /(</>i, Tg, Ty) gives the change of Ty due to 
the coupling of the S-beats to the V-beats . Each 

of the 2 or 3 S-beats that may appear between consec- 
utive (expressed or blocked) V-beats iteratively changes 
the position of the next V-beat depending on the timing 
of the S-beats as shown in Fig. |4|b. These changes add 
up to T s ,Ty) ra. 
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(a) Deterministic Model 



(b) Data 



(c) Model with Noise+Coupling 
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FIG. 3. Histograms of the number of intervening beats (NIB) on a logarithmic scale for Ts = 0.51 s, 0.55 s, and 0.61 s, from 
top to bottom. For Ty =1.67 s, the corresponding ratios of the periods are Ty /Ts = 3.27, 3.03, and 2.73. In the middle panel 
the count for NIB=2 is 2465 in (a), 677 in (b), and 686 in (c). (a) Simulation of the deterministic model with constant Ts, 
Ty = 1-67 s. The average S-beat interval Ts varies from 0.47 s to 0.64 s based on the dataset we study. For each histogram, the 
number of S-beats corresponds to the number of S-beats in the data with the same average Ts- (b) NIB- histogram for the data. 
The value for Ts at each V-beat is the average S-beat interval in a window of 20 beats centered at the V-beat. All V-beats in 
this window are replaced by the corresponding blocked S-beats which are placed in the middle of the two neighboring S-beats. 
The histogram of the NIB-values is computed from all sequences of the record that have the same mean Ts. For Ts = 0.55 s 
(TV/Ts=3.03), we find NIB-values 2, 5, 8, 11, . . . , 3n — 1 that do not appear in the deterministic model, (c) Simulation with 
coupling (Fig. Qs) and noise. Instead of a constant Ts, we use the S-beat intervals of the data, and we add to Ty Gaussian 
white noise with a standard deviation a = 0.07 s. Note that the intrinsic period Ty = 1.76 s is larger than the apparent period 
Ty ~ 1.67s because of the coupling. 



In order to reproduce Fig. we have chosen a coupling 
that leads, for T v /T s « 3, to a fixed point in Eq. (§), 
such that the V-beats always fall in the refractory period 
of the S-beat and thus are always blocked. However, the 
noise term leads to a dispersion of the phases of the V- 
oscillator such that some of the V-beats are expressed. 
Thus, both the coupling and the stochastic term interact 
to generate the dynamics. The effects of the coupling 
and the noise for Ty /Ts ^ 3 are represented schemati- 
cally in Fig. ^b. Fig. |3|c shows how the distributions of 
the NIB-values change when both coupling and noise are 
included in the simulation. The model now reproduces 
quantitatively the data. For T s ~ 0.55 s (T v /T s w 3) the 
number of the V-beats is reduced by the coupling, and 
the NIB-sequence 2,5,8,... is generated by the noise. As 
a consequence of the random noise moving the V-beats 
randomly in and out of the refractory time, the blocking 
mechanism gives rise to a discrete Poisson process lead- 
ing to an approximately exponential fall-off of the peak 
heights of the occurrences of the NIB values 2,5,8 
Finally, the model gives an accurate estimate of the frac- 
tion of V-beats as a function of T$ shown in Fig. ||a. 

In order to estimate the magnitude of the noise, we 
calculate the cross-correlation between the numbers of 
occurrences of each NIB value in the model and the clin- 
ical data pq] . Fig. [f| shows the correlation as a function 
of the standard deviation a of the noise. The correlation 
function has a maximum at a *=s 0.07 s. This value of 
a also best reproduces the broadness of the peaks in the 
distribution of the interventricular time intervals (Fig. ^) . 
Further, the model reproduces the asymmetrical form of 
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FIG. 4. (a) The fraction ny/N plotted against the ratio 
Ty/Ts- The theoretical curve given by Eq. (0), with 8 as in 
Fig. |^a, is reproduced by the model without coupling, and re- 
mains unchanged when noise is added. The data deviate from 
this curve at the ratio Ty/Ts = 3, where we find less than half 
of the predicted V-beats. The model with coupling and noise 
reproduces this behavior, (b) The coupling between the two 
oscillators. The change ATV of the period of the V-oscillator 
as a function of the ratio of the time tvs between the last 
V-beat and an S-beat, and Ty. The coupling shortens the 
intrinsic Ty to the apparent Ty. 
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the peaks in the histogram (inset of Fig. [T]). In the sim- 
ulation, the coupling splits the peaks into two subpeaks 
giving them an asymmetrical appearance (not shown). 

In this work we analyzed the patterns of abnormal 
heartbeats in a 12 hour record from a single patient and 
proposed that the dynamics results from a combination 
of deterministic and stochastic mechanisms. Quantita- 
tive comparison between predictions of the model and 
the clinical data shows best agreement for an optimal 
level of noise in the model. 

Our approach is in contrast to standard approaches |l]] 
in which crude measures, such as the average numbers of 
abnormal heart beats per unit time, are used for clini- 
cal assessment. We believe that distinctive dynamics are 
associated with different mechanisms, and hence differ- 
ent therapeutic strategies. Thus, the detailed program 
of analysis applied to one clinical record in this Letter is 
essential to better classify cardiac arrhythmias based on 
the underlying mechanisms and resulting dynamics. 
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FIG. 5. The cross-correlation between data and model as 
a function of noise. We cross-correlate the histograms in 
Fig. HbjC for all values of Ts from 0.48-0.62 s in steps of 
0.01s ||. 
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